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1 INTRODUCTION 


ABSTRACT 

The VST Photometric Ha Survey of the Southern Galactic Plane and Bulge (VPHAS+) is 
surveying the southern Milky Way in u, g, r, i and Hæ at ~1 arcsec angular resolution. Its 
footprint spans the Galactic latitude range —5° < b < +5° at all longitudes south of the 
celestial equator. Extensions around the Galactic Centre to Galactic latitudes +10° bring in 
much of the Galactic bulge. This European Southern Observatory public survey, begun on 
2011 December 28, reaches down to ~20th magnitude (100) and will provide single-epoch 
digital optical photometry for ~300 million stars. The observing strategy and data pipelining 
are described, and an appraisal of the segmented narrow-band Ha filter in use is presented. 
Using model atmospheres and library spectra, we compute main-sequence (u — g), (g — r), 
(r — i) and (r — Haq) stellar colours in the Vega system. We report on a preliminary validation 
of the photometry using test data obtained from two pointings overlapping the Sloan Digital 
Sky Survey. An example of the (u — g, g — r) and (r — Ha, r — i) diagrams for a full VPHAS+ 
survey field is given. Attention is drawn to the opportunities for studies of compact nebulae 
and nebular morphologies that arise from the image quality being achieved. The value of the u 
band as the means to identify planetary-nebula central stars is demonstrated by the discovery 
of the central star of NGC 2899 in survey data. Thanks to its excellent imaging performance, 
the VLT Survey Telescope (VST)/OmegaCam combination used by this survey is a perfect 
vehicle for automated searches for reddened early-type stars, and will allow the discovery and 
analysis of compact binaries, white dwarfs and transient sources. 


Key words: surveys—stars: general—Hu regions—planetary nebulae: general- Galaxy: 
bulge — Galaxy: disc. 


way of developing and testing models for these crucial evolutionary 
stages. 


The Ha emission line is well known as a tracer of diffuse ion- 
ized nebulae and as a marker of pre- or post-main-sequence status 
among spatially unresolved stellar sources. Since these objects — 
both nebulae and stars — represent relatively short-lived phases of 
evolution, they amount to a minority population in a mature galaxy 
like our own. Their relative scarcity has in the past stood in the 
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In the southern hemisphere, the search for planetary nebulae 
(PNe) has been served well by Ha imaging surveys carried out by 
the UK Schmidt Telescope (Parker et al. 2005, 2006 and other more 
recent works). Nevertheless, the VST Photometric Ha Survey of the 
Southern Galactic Plane and Bulge (VPHAS-+-) will have a decisive 
impact on studies of complex or smaller nebulae of all types, ranging 
from optically detectable ultra-compact and compact H 11 regions, 
to nebulae around young stellar objects (including associated jets 
and Herbig-Haro objects), through PNe, to extended emission from 
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D-type symbiotic stars and supernova remnants. The superb spatial 
resolution, dynamic range and likely photometric accuracy of the 
VPHAS-+ images warrant a step forward in our knowledge of the 
population and detailed characteristics of these object classes. 

For southern point sources with line emission, the situation is very 
different: there has been little updating of the available catalogues 
since the work of Stephenson & Sanduleak (1971) that was limited 
to a depth of 12th magnitude. The major groups of emission line 
stars that remain as challenges to our understanding include all 
types of massive star [O stars, supergiants, luminous blue variables 
(LBV), Wolf-Rayet (WR) stars, various types of Be star], post- 
asymptotic giant branch stars, pre-main-sequence stars at all masses, 
active stars and compact interacting binaries. Within the disc of the 
Milky Way, the available samples of these objects are typically 
modest and heterogeneous. Fixing this deficit via a uniform search 
of the Galactic plane for these rare object classes motivated the 
photometric Ha survey of the southern Galactic plane, first proposed 
for the VLT Survey Telescope (VST) in 2004. This paper describes 
the realization of this European Southern Observatory (ESO) public 
survey, now known as the VST Photometric Ha Survey of the 
Southern Galactic Plane and Bulge (VPHAS-+). 

When first proposed, VPHAS (without the plus sign) was envis- 
aged as the counterpart to the INT/WFC Photometric Ha Survey 
of the Northern Galactic Plane (IPHAS; Drew et al. 2005), which 
had begun in 2003 August. IPHAS is a digital imaging survey made 
up of Sloan r, i and narrow-band Hø exposures, reaching to ~20th 
magnitude, which takes in all Galactic longitudes north of the ce- 
lestial equator in the latitude range —5° < b < +5°. This is all 
but complete with the new release of a catalogue of ~200 million 
unique objects drawn from 92 percent of the survey’s footprint 
(Barentsen et al., in preparation). During the initial VST public 
survey review process, it was agreed that VPHAS should broaden 
in scope to also incorporate the Sloan u and g bands (proposed 
for a separate survey), which are particularly useful in picking out 
observation block (OB) stars, white dwarfs and other blue-excess 
objects. With this upgrade to five bands, the renamed VPHAS-+ be- 
came an all-purpose digital optical survey of the southern Galactic 
plane, capable of delivering data at a spatial resolution of ~1 arcsec 
or better. As well as fulfilling the role of southern counterpart of 
IPHAS, VPHAS-+ is also the counterpart to UVEX, the UV-Excess 
Survey of the Northern Galactic Plane (Groot et al. 2009) that, at 
the time of writing, continues on the Isaac Newton Telescope in La 
Palma. 

The final augmentation of the VPHAS-+ survey footprint came in 
2010 on expanding its footprint to match that of the similarly high- 
spatial-resolution near-infrared (NIR) survey, VISTA Variables in 
the Via Lactea (VVV; Minniti et al. 2010). The survey footprint 
now includes the Galactic bulge to a latitude of |b| < 10°, across the 
longitude range — 10° < £ < +10°. The VVV z, Y, J, H and K; survey 
of much of the bulge and inner Galactic disc is already complete. 

VPHAS+ is poised to become the homogeneous digital optical 
imaging survey of the Galactic plane and bulge, at ~1 arcsec angular 
resolution, that will provide a uniform database of stellar spectral en- 
ergy distributions (SEDs), from which a range of colour-magnitude 
and colour—colour diagrams of well-established utility can be de- 
rived. Once calibrated to the expected precision of 2-3 per cent, like 
its northern counterparts, VPHAS-+ will be quantitatively far supe- 
rior to the photographic surveys of the last century, and will offer 
significant added value in the form of calibrated narrow-band Ha 
data. Some of the science enabled is illustrated by the studies that 
the northern surveys, IPHAS and UVEX, have already stimulated. 
These have included a number of works reporting the discovery of 
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emission line stars, ranging from young objects (e.g. Vink et al. 
2008; Valdivieso et al. 2009; Barentsen et al. 2011; Raddi et al. 
2013) to evolved object classes such as symbiotic and cataclysmic 
binaries and compact PNe (e.g. Witham et al. 2006; Wesson et al. 
2008; Viironen et al. 2009a; Corradi et al. 2010). The diagnostic 
power to be expected from the blue u and g bands has been ap- 
praised in a series of papers presenting UVEX and early follow-up 
spectroscopy (Groot et al. 2009; Verbeek et al. 2012a,b). 

Any comprehensive survey of the Galactic plane clearly targets 
the main mass component of our own Galaxy, made up of stars, 
gas and dust. When a spatial resolution of 1 arcsec is combined 
with wide area coverage spanning hundreds of square degrees, as 
in IPHAS, UVEX and VPHAS-+, it becomes possible to exploit the 
stellar photometry adaptively to solve for the distribution of both 
the stars and dust making up the optically accessible Galactic disc 
— by means that are similar to those already attempted, based on 
the NIR 2MASS survey (Drimmel & Spergel 2001; Marshall et al. 
2006). Methods to achieve three-dimensional mapping of this kind, 
now incorporating the direct sensitivity of the Ha narrow band to 
stellar intrinsic colour, are starting to take shape (Sale et al. 2009; 
Sale 2012). This development coincides with the approach of the 
operations phase of Europe’s next major astrometric mission, Gaia. 
Indeed, the rich dynamical picture that Gaia will build of the Milky 
Way over the next decade will be very effectively complemented 
by stellar energy distributions measured for millions of stars from 
the current generation of ground-based optical, and NIR, wide-field 
surveys. VPHAS+, with its haul of photometry in five optical bands 
on ~300 million objects, is set to take its place as one of them. 

This paper presents the main features of VPHAS+, including 
a description of its execution, the data processing and the nature 
of the photometric colour information it provides. We begin in the 
next section with a presentation of the observing strategy and the 
data reduction techniques in use. We then turn to a description and 
evaluation of the narrow-band Ha filter procured for this survey, in 
Section 3. Following this, in Section 4, we present tailored synthetic 
photometry of main-sequence and giant stars that provides insights 
into the photometric diagrams that may be generated from survey 
data. An exercise in photometric validation is described in Section 5 
in which Sloan Digital Sky Survey (SDSS) data are compared with 
VST observations. The scene is then set for an example of VPHAS+ 
photometry extracted across the entire square-degree footprint of a 
single survey field (Section 6). In Sections 7 and 8, we outline the 
applications of VPHAS-+ to spatially resolved nebular astrophysics 
and in the time domain. The paper ends in Section 9 with a summary, 
examples of early survey exploitation and a forward look to the first 
major data release. 


2 SURVEY OBSERVATIONS AND DATA 
PROCESSING 


2.1 VPHAS+ specification 


The footprint of the survey is shown in Fig. 1. The OmegaCAM 
imager (Kuijken 2011) on the VST provides a field size of a full 
square degree, captured on a 4 x 8 CCD mosaic. After allowing 
for some modest overlap between adjacent fields, we arrived at a 
set of 2269 field centres that will cover the desired Galactic latitude 
band —5° < b < +5° at all southern hemisphere Galactic longitudes, 
as well as incorporate the Galactic bulge extensions to —10° < 
b < +10° near the Galactic Centre. The survey footprint extends 
across the celestial equator by a degree or two to achieve an overlap 
with the northern hemisphere surveys IPHAS and UVEX of ~100 
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Figure 1. The VPHAS-+ survey footprint plotted in Galactic coordinates. All 2269 fields are shown in outline. Different colours, as specified in the key, are 
used to identify the observations obtained for each field by 2014 January | — essentially 2 years after the start of data taking. 


square degrees altogether. This is to create the opportunity for some 
direct photometric cross-calibration. 

The target depth of the survey is to reach to at least ~20th mag- 
nitude, at 100, in each of the Sloan u, g, r and i broad-band filters 
and narrow-band Ha. The bright limit consistent with this goal is 
typically 12—13th magnitude. Presently, all VWPHAS+ photometric 
magnitudes are expressed in the Vega system. The original con- 
cept was to collect the data in all five bands contemporaneously, 
in order to build a uniform library of snapshot photometric SEDs 
for 200 million or more stars. Practical constraints have modified 
this to the extent that the blue filters (u, g) are observed separately 
from the reddest (i and Ha), with the r band serving as a linking 
reference that is observed both with u, g and i, Ha. The aim is also 
to keep the spatial resolution close to 1 arcsec. OmegaCAM and 
the Paranal site are well suited to this in that the camera pixel size 
is 0.21 arcsec, projected on sky, and the median seeing achieved is 
better than 1 arcsec (on occasion falling to as little as 0.6 arcsec). 

As a means to obtaining better quality control, and to ensure that 
only a minimal fraction of the survey footprint is missed due to the 
pattern of gaps between the CCDs in the camera mosaic, every field 
is imaged at two or three offset pointings. This strategy has been 
carried over from IPHAS and UVEX in the northern hemisphere, 
and has the consequence that the majority of imaged objects will be 
detected twice some minutes apart. In the r band, there will be two 
arbitrarily separated epochs of data, with typically two detections 
at each of them (i.e. four altogether). 


2.2 VPHAS-+ observations 


The VST is a service-observing facility, with all programmes 
queued for execution as and when the ambient conditions meet 
programme requirements. VPHAS-+ survey field acquisition began 
on 2011 December 28. Normally the constraint set includes a seeing 
upper bound of 1.2 arcsec: this is only set at a lower, more stringent 
value for fields expected to present a particularly high density of 
sources (e.g. in the southern bulge). In order that the seeing achieved 


MNRAS 440, 2036-2058 (2014) 


in the u band is not greatly different from that in 7 at the opposite 
end of the optical range, it is advantageous to separate acquisition of 
blue data from red — hence a split between ‘blue’ (u, g, 7) and ‘red’ 
(Ha, r, i) observing blocks has been implemented. This split also 
permits the use of different moon distance and phase constraints, 
such that blue data are obtained when the moon is less than half full 
at an angular separation of not less than 60°, while the limits for 
red data are set at 0.7 moon illumination and a minimum angle of 
50°. Avoiding bright-moon conditions is important in order to limit 
the amount of moonlight mixed in with diffuse Ha emission in the 
reduced images. No requirement has been placed on the time elaps- 
ing between acquisition of blue and red data. However, the more 
forgiving constraints on the acquisition of the latter have meant that 
these are typically executed sooner than the former, with the result 
that many more fields have red data already than blue data (see 
Fig. 1). In all cases, the final constraint is that the sky is required to 
be clear, if not necessarily fully photometric. 

An impressive feature of the camera, OmegaCAM, is its potential 
to deliver remarkably undistorted point-source images all the way 
across the 1° field of view. To realize this, it is critical that the 
VST has an actively controlled primary. The operational price for 
this, at the present time, is that image analysis and correction have 
to be carried out at every filter change or after longer slews. The 
overhead added by this is about 3 min. To reduce the impact of 
this, observations of sets of three neighbouring fields are scheduled 
together, so that the image analysis needs to only take place every 
15-30 min — not much more often than would be essential, in any 
case, to compensate for the telescope’s tracking movement. As a 
result, ‘contemporaneous’ in the context of VPHAS+ data taking 
means that all three blue, or red, filters are typically exposed within 
40-50 min of each other (cf. IPHAS, where the more compact 
camera allows much faster operation, bringing this elapsed time 
down to under 10 min). However, the time difference between the 
blue and red observing blocks for a given field, i.e. the u/g/r data 
collection and Ha/r/i data collection, can be anything from a few 
hours to more than a year. 
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Figure 2. The transmission profiles of the Sloan u, g, r and i, and narrow- 
band Ha filters used in all VPHAS+ observations. The r and Hæ profiles 
are shown as a dashed line and in red, respectively, just to clearly distinguish 
them from each other and the i band. Each profile has been multiplied by 
the CCD response function and a model of atmospheric throughput (Patat 
et al. 2011). The very nearly grey losses due to the telescope optics (a further 
scaling of approximately 0.6) have not been folded in. 


Table 1. Observations obtained per survey field. The 
median seeing quoted is derived from data in hand by 
2013 December. 


Filter Exposure No. of Median seeing 
time (s) offsets (arcsec) 
Blue observation blocks 
u 150 2 1.01 
g 40° 3 0.88 
r 25 2 0.80 
Red observation blocks 
Ha 120 3 0.84 
r 25 2 0.82 
i 25 2 0.77 


“Up to 2013 February 19, g exposure times were 30 s. 


We provide a reminder of the passbands of the Sloan broad-band 
filters in use, along with that of the narrow-band Ha, in Fig. 2. They 
are shown scaled by a typical CCD response function and a model 
of the atmospheric transmission (Patat et al. 2011). The exposure 
times used for the different filters, the number of exposures in each 
and the median seeing achieved, up to 2013 December, are set out 
in Table 1. From early on, during the commissioning phase of the 
telescope, it became clear that tracking is usually good enough that 
even the 150 s u exposures, our longest, do not have to be guided in 
most circumstances. Indeed, experience is showing that it is safer 
to rely on the tracking, rather than the autoguider, to maintain good 
image quality in the most dense star fields. 

The pattern of offsets used for each field is illustrated in Fig. 3. 
The shifts are relatively large, with the outer pointings differing 
by —588 arcsec in the RA direction and +660 arcsec in declination. 
The choices made have largely been driven by the characteristics 
of the narrow-band Ha filter (discussed below in Section 3), but 
they also convey the advantage of greatly increasing the overlaps 
between neighbouring fields. Just the two outermost pointings are 
used when exposing the u, r and i filters. This leaves 0.4 per cent of 
the survey footprint unexposed. This changes to complete coverage 
on including the third intermediate offset, as is the policy for the 
Ha and g filters. 
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Figure 3. An illustration of the VPHAS+ offset pattern as it applies to 
the segmented Ha filter with extra vignetting due to the T-bars separating 
the four segments. The first pointing is to lower right, as drawn — a con- 
servative estimate of the exposed unvignetted area is shown in black. The 
exposed/unvignetted areas of the second and third pointings are shown in 
green and blue, respectively. The vertical and horizontal scales are numbered 
in pixels (RA increasing to the right, declination increasing upwards). 


In accordance with ESO’s standard procedures, data are evalu- 
ated soon after collection by Paranal staff and graded before transfer 
to the archive in Garching and to the Cambridge Astronomy Survey 
Unit (CASU) in Cambridge. If the applied constraints are signifi- 
cantly violated, the observation block is returned to the queue. 


2.3 Data pipeline 
2.3.1 Initial processing 


From 2012 February 29, raw VST data have been routinely trans- 
ferred from Paranal to Garching over the internet. For each ob- 
servation, the imaging data are stored in a Multi-Extension FITS 
file with a primary header describing the overall characteristics of 
the observation (pointing, filter, exposure time, etc.) and 32 image 
extensions, corresponding to each of the CCD detectors, with fur- 
ther detector-level information in the secondary headers. The 32-bit 
integer raw data files are Rice-compressed at source using loss- 
less compression (e.g. Sabbey, Coppi & Oemler 1998). The files 
are then checked and ingested into the ESO raw data archive in 
Garching. As soon as the data for any given night become avail- 
able, they are automatically transferred to the CASU for further 
checks and subsequent processing. The VST web pages at CASU 
provide an external interface for both monitoring processing sta- 
tus (http://casu.ast.cam.ac.uk/versus/data-processing/) and overall 
survey progress and access (http://casu.ast.cam.ac.uk/vstsp/). 

The processing sequence is similar to that used for the IPHAS 
survey of the northern Galactic plane (e.g. Gonzalez-Solares et al. 
2008), while the higher level control software is based on that devel- 
oped for the VISTA Data Flow System (VDFS; Irwin et al. 2004). 
Here we briefly outline the processing steps illustrated in Fig. 4, 
emphasizing the main differences relative to the current VDFS stan- 
dard. A more detailed description of the VST processing pipeline 
is currently in preparation (Kupcu Yoldas et al., in preparation). 

Science images are first debiased. Full two-dimensional bias re- 
moval is necessary due to amplifier glow during readout being 
present in some detectors. The master bias frames are updated daily 
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Figure 4. Flowchart identifying the main VST data processing steps. At the present time, all steps up to the reduced single-band catalogue are undertaken at 


CASU. Band merging is performed at the University of Hertfordshire. 


from calibration files taken as part of the operational cycle. The 
OmegaCAM detectors are linear to better than 1 per cent over their 
usable dynamic range removing the need for a linearity correction. 
Hence, this stage in the pipeline processing (Fig. 4), although part 
of the pipeline architecture, is currently bypassed. 

Flat-field images in each band are constructed by combining a 
series of twilight sky flats obtained in bright sky conditions. The 
time-scale to obtain sequences of these for all deployed filters is 
typically one to two weeks. So as to adequately trace the variations 
in the pattern and level of scattered light in these flats, the master 
flats derived from them are updated on a monthly cycle (how these 
are corrected for scattered light is described in Section 2.3.3). Four 
of the detectors, in extensions 29—32, suffer from inter-detector 
cross-talk, whereby saturated bright stars in one detector can cause 
noticeable positive or negative low-level (~0.1 percent) ghost ob- 
jects in adjacent detectors. The impact of these is minimized in the 
pipeline by applying a pre-tabulated cross-talk correction matrix to 
each of the affected images. 

The flat-field sequences plus bad pixel masks are used to generate 
the confidence (weight) maps (e.g. Irwin et al. 2004) used later 
during catalogue generation and any subsequent image stacking or 
large-area mosaicking. After flat-fielding, science images generally 
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have well-behaved sky backgrounds which makes subsequent image 
processing straightforward. Where direct scattered light is present in 
them, it is an additive phenomenon that is dealt with automatically 
during object catalogue generation. 

The redder passbands used in VST observations, in this case the 
i band, show fringing patterns at ~2 per cent of the sky background 
level. Defringing is done using a standard CASU procedure (Irwin 
& Lewis 2001). The fringe frames used are derived from other 
VST public survey data taken as close as possible in time at higher 
Galactic latitude. This approach works because the fringe pattern 
induced by sky emission lines at Paranal is quite stable over long 
periods. The fringe frames are automatically scaled and subtracted 
from each science image reducing the residual fringing level to well 
below the sky noise. 

Catalogue generation is based on mcore! (Irwin 1985) and makes 
direct use of the confidence maps, derived from the flat-fields, to 
suitably weight down unreliable parts of the images. This step in- 
cludes object detection, parameterization and morphological clas- 
sification, together with generation of a range of quality-control 


' Software publicly available from http://casu.ast.cam.ac.uk. 
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Figure 5. The left-hand panel shows an example of the deduced scattered light component present in 2013 June r-band data. The right-hand panel shows 
the corrected outcome. These maps were constructed using 286 318 APASS object matches over the square-degree field. Before correction, scattered light 
gradients amounting to a 20-25 per cent variation from corners to centre are present. This flattens to around +2 per cent. 


information. Because of the extensive presence of diffuse emission 
throughout the southern Galactic plane, particularly in Ha, a version 
of each affected image is cleaned of nebulosity using the NEBULISER” 
for the purpose of catalogue generation only. This achieves a more 
careful removal of background and ultimately leads to more com- 
plete and, on average, more faithful object detection than in the 
absence of this step. 

With object catalogues available for every VPHAS+ survey im- 
age, it is then possible to improve the rough World Coordinate 
System (WCS) based on the telescope pointing and general system 
characteristics. The WCS is progressively refined using matches be- 
tween detected objects and the 2MASS catalogue (Skrutskie et al. 
2006). Despite the large field of view, the VST focal plane is almost 
free of distortion, and a standard tangent plane projection yields 
residual systematics of ~25 mas over the entire field. 


2.3.2 Photometric calibration 


Provisional photometric calibration is based on a series of standard 
star fields observed each night (e.g. Landolt 1992). For each night, a 
zero-point and error estimate using the observations of all the stan- 
dard fields in each filter is derived. The flat-fielding stage nominally 
places all detectors on a common internal gain system implying, in 
principle, that a single zero-point suffices to characterize the whole 
focal surface. Colour equations are used to transform between the 
passbands in use on the VST and the Johnson—Cousins system of 
the published standard star photometry. The calibration is currently 
in a VST system that uses the SED of Vega as the zero-colour, 
almost zero-magnitude, reference object. 

The u-band data are the most challenging to calibrate. As this 
part of Vega’s spectrum, and also the average standard star plus the 
detector response, is falling rapidly, it would be surprising if there 
were no offsets in u due to non-linearities in the required colour 
transforms and, perhaps, to degenerate colour transforms for hotter 
stars. Early experience of working with u data does indeed suggest 


? Software publicly available from http://casu.ast.cam.ac.uk, see also Irwin 
(2010). 


that offsets of up to a few tenths of a magnitude are sometimes 
present (see Sections 5 and 6). 

The colour transforms currently in use to define the VPHAS 
internal system are given below: 


uyst = U + 0.035 (U — B) 
gvst = B — 0.405 (B — V) 
rvst = R + 0.255 (V — R) 
iysr = I + 0.215 (R — I) 
Haysr = R + 0.025 (V — R). 


The transform for the narrow-band Ha is an approximate initial 
solution needed for the subsequent illumination correction stage. 
At catalogue band merging, this is superseded (see Section 4). 


2.3.3 Illumination correction 


The main difficulty in deriving an accurate photometric calibration 
over the 1° field arises from the multiplicative systematics caused 
by scattered light in the flat-fields. The VST (at least up to the 
introduction of baffles early in 2014) has proved particularly sus- 
ceptible to variable scattered light. Its impact has varied from month 
to month depending on conditions prevalent at the time the flat-field 
sequences were taken. An illustration of the amount and character 
of the master-flat correction required is provided in Fig. 5. 

The scattered light is made up of multiple components with 
different symmetries and scales. These range from ~10 arcsec 
with x-y rectangular symmetry, e.g. due to scattering off mask- 
ing strips above the CCD readout edges, to large fractions of the 
field due to radial concentration of light in the optics and to non- 
astronomical scattered light entering obliquely in flat-field frames. 
After some experimentation, and external verification, we found that 
the AAVSO Photometric All-Sky Survey (APASS) g,r,i catalogues 
(http://www.aavso.org/apass) provide a reliable working solution 
to the illumination correction problem inherent in VST data (see 
Fig. 5). These catalogues also provide an independent overall pho- 
tometric calibration tied to the SDSS AB magnitude system and 
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Figure 6. Photometric errors in VPHAS+ data as a function of magnitude. All data are drawn from an area of ~0.2 square degrees in field 1679, 
positioned ~20 arcmin E of Westerlund 2. The left-hand panel refers to blue OB data, while the right refers to red OB data. The coloured histogram in each 
component plot shows the mean absolute deviation of the magnitude difference, (mı — mz), per 0.25 mag bin, while the bin means of the suitably corrected 
pipeline estimate of the random error on this difference are in black. The 100 and 50 magnitude limits are specified in the parentheses next to each filter name. 
The rising observed mean deviation seen in g at the bright end is due to the onset of saturation. 


will be used in future updates to define an alternative finer grained 
temporal AB magnitude zero-point. 

All filters used are treated in the same manner with colour equa- 
tions set up to define transformations between the APASS g,r,i 
SDSS-like calibration and the VPHAS-+ u,g,r,i,Hæ internal system. 

Illumination corrections are re-derived for each filter once a 
month. Application of these corrections via the master flats re- 
duces the residual systematics across the entire field to below the 
1 percent level for the broad-band filters and to within 2 percent 
for the segmented narrow-band Ha, except in vignetted regions (see 
Section 3). 


2.3.4 Quality control 


In addition to the usual VDFS quality-control monitoring of aver- 
age stellar seeing and ellipticity, sky surface brightness and noise 
properties, we have also initiated a more detailed analysis of the 
image properties based on inter-detector comparisons. 

The well-aligned coplanar detector array coupled with the curved 
focal surface is extremely sensitive to imperfections in focus which 
are relatively easy to detect using the detector-level average seeing 
measurement variation available for each of the 32 detectors. Like- 
wise, the variation in average stellar ellipticity from each detector 
over the field is used to monitor rotator angle tracking problems. 

All of this information can be used in addition to the OB grades 
provided by ESO and is incorporated within all data product files 
and also the progress data base. 
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2.4 Limiting magnitudes and errors 


The present convention for VPHAS+ and this paper is that all 
magnitudes are expressed in the Vega system, which imposes zero 
intrinsic colour for AO stars. The 5ø limiting magnitudes commonly 
achieved per exposure range from 20.5—21.0 for Ha up to 22.2—22.7 
for Sloan g. The 100 limits are about 1 mag brighter. 

Every source flux or magnitude determined via the pipeline has 
a formal error associated with it. We provide an example of how 
these compare with empirical magnitude differences, by extract- 
ing a sample of stars from a 0.25-square-degree catalogue, cut out 
from the survey field including Westerlund 2 (field 1679, see also 
Section 6) in order to examine the pattern of errors (Fig. 6). The sky 
area chosen is offset from the cluster to the east by ~20 arcmin and 
exhibits moderate diffuse ionized nebulosity. In the southern plane, 
the presence of some nebulosity, particularly affecting r and Ha 
exposures, is more the rule than the exception. Sources classified 
as probable stars in both g and r (blue-filter set, left-hand panel) 
and r and i (red-filter set, right-hand panel) in two consecutive off- 
set exposures have been selected. The selection also required that 
each extracted magnitude was unaffected by vignetting and bad pix- 
els (confidence level >95). This step is particularly important for 
Ha given the extra vignetting introduced by the cross bars of the 
segmented filter (see Section 3). 

The faintest stars that might have been included in the plots for i 
and blue r (top row in Fig. 6) are absent because of a requirement 
that every included source should also be picked up in, respectively, 
red r and g. Fainter objects than the apparent limits certainly exist 
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in these bands. This feature follows directly from the typically red 
colours of Galactic plane stars at magnitudes fainter than ~13 that 
are the target of this survey. For the same reason, it is not uncommon 
for the u-band source counts to be one or more orders of magnitude 
lower than those of the į band. The role of the u band is to pick out 
the unusual rather than to characterize the routine. 

To bring out the systematic effects present, the specific com- 
parison made in Fig. 6 is between the bin means of the absolute 
magnitude differences, |n; — m — ô|, between the two exposures 
and the expected random error on the difference derived from the 
pipeline rms errors on the individual magnitude measurements. The 
quantity, 5, is the median magnitude difference computed from all 
bright stars down to 18th magnitude (r, i and Ha) or 19th magnitude 
(u and g). This was small in all cases — the largest value being 0.011 
for u. On the other hand, the correction applied to the pipeline errors 
was, first, to multiply the single-measurement magnitude error by 
/2 to give the rms error on the (mı — m2) difference, and then 
to multiply by ./2/7« in order to convert the measure of dispersion 
from rms to a mean deviation. 

At magnitudes brighter than 18-19 in Fig. 6, the scatter in the 
empirical results can be seen to be appreciably greater than that 
‘predicted’ for the random component by the pipeline. The scale of 
the difference indicates that a further error component of ~0.01- 
0.02 mag is present. The amount and filter dependence of these 
levels of error are entirely consistent with the uncertainties esti- 
mated above for the flat-field and illumination corrections: as noted 
in Section 2.3.3, the VST is presently prone to quite high and vari- 
able levels of scattered light. Practical remedies for this are under 
consideration by ESO — when implemented these should tighten up 
the error budget. 

The enhanced mean magnitude difference seen for g < 13 in 
Fig. 6 is typical of what is seen as saturation effects begin to set 
in. For all the other bands, in this example, saturation sets in at 
magnitudes a little brighter than 12. A safe working assumption 
across VPHAS+ would be that saturation is never troublesome at 
magnitudes fainter than 13, but always an issue for magnitudes 
brighter than 12. 

Fig. 6 identifies the 5o and 100 magnitude limits for each filter 
achieved in this representative example. The seeing at the time of 
these observations, as measured in the pipeline, ranged from 0.8 to 
1 arcsec (cf. Table 1). 


3 THE NARROW-BAND Ha FILTER 


3.1 Overview 


A filter required to select a narrow band across a large 27x27 cm? 
image plane is a challenging fabrication problem. At the telescope, 
the filter in use for VPHAS-+ is known as NB-659. At the time it was 
commissioned in 2006, the purchase of a single-piece narrow-band 
filter was offered only by one supplier and was well beyond budget. 
This left the four-segment option as the achievable alternative. 
The Ha filter was constructed based on a specification supplied 
by the OmegaCAM consortium, setting as goal a central wavelength 
(CWL) of 6588 A and a bandpass of 107 A. It was delivered in the 
summer of 2009, and was shortly thereafter tested at the University 
of Munich Observatory, using the optical lab set up by the Omega- 
CAM consortium for filter testing. A photo of the filter at that time 
is shown in Fig. 7. The transmission of each filter segment was 
measured at 21 positions forming a coarse radial pattern (Fig. 8) 
using a monochromator beam adjusted to emulate the f-ratio 5.5 
VST/OmegaCAM optical system. The logic of the chosen measure- 
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Figure 7. The segmented Ha filter, photographed in the lab soon after 
receipt and just prior to measuring its transmission. The filters to either side 


of the 2x2 array of Ha segments, transmit as the r band and cover the guide 
CCDs. 
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Figure 8. A map of the positions within each filter segment at which the 
transmission was measured. The colour codings are used again in later 
figures to distinguish the central positions (black, crosses, sampling roughly 
30 per cent of the segment area), intermediate radii (blue, encircled crosses, 
sampling ~half the area) and corners (red crosses, 15 per cent of the area — 
of which almost half is lost to vignetting). The dashed lines define the limits 
of the strips 2 arcmin wide that experience any vignetting due to the filter 
T-bars. They are drawn here as for segment D. 


ment pattern is to give a good sampling of the dominantly radial 
variation of the transmission profile due to the turntable rotation 
in the filter coating chamber. The diameter of the monochroma- 
tor beam used in the measurements was 4-5 mm. This is a more 
compact beam than that of starlight at the telescope, which fills a 
spotsize of up to 12 mm on passing through the filter out of focus. 
Consequently, the actual performance will be a somewhat a really 
smoothed version of the performance revealed by the lab measure- 
ments and their subsequent simulation. The filter was shipped to 
Paranal and VST in the spring of 2011, after some final selective 
remeasuring. These confirmed that there had been no discernible 
bandpass changes in store since delivery almost 2 years earlier. 

At the time the monochromator measurements were made, a 
segment naming scheme was put in place (segments A, B, C and 
D) which is re-used here. Presently, the filter is housed in magazine 
B of OmegaCAM, which means that in terms of the view of the 
sky, segment A spans the SW section of the image plane, B the SE, 
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Table 2. Summary of filter segment properties. 


Segment Sky quadrant CCDs covered Centre, mean, corner CWL Mean integrated 
(A) throughput (A) 
A SW 1-8 6580.2, 6585.4, 6595.3 98.64 
B SE 17-24 6596.1, 6585.4, 6578.9 103.38 
C NE 25-32 6582.8, 6591.9, 6599.8 99.74 
D NW 9-16 6581.7, 6594.3, 6603.8 99.49 
1.0 T Table 3. Mean transmission for NB-659. 
0.9 - = 
Wavelength Transmission Wavelength Transmission 
a, ] (A) (A) 
0.7 4 
a 6456.3 0.000 6591.5 0.962 
p eF 7 6461.3 0.001 6696.5 0.961 
g 0.5 H 4 6465.8 0.001 6601.0 0.960 
a zal | 6470.8 0.002 6616.0 0.955 
x i 6475.9 0.002 6611.0 0.945 
ig a 6480.9 0.003 6616.1 0.928 
aale ] 6485.9 0.005 6621.1 0.896 
6490.9 0.008 6626.1 0.839 
Saki J 6496.0 0.012 6631.1 0.736 
n f f fi f fi i L L 6501.0 0.020 6636.2 0.609 
6460 6500 6540 6580 6620 6660 6700 6740 6506.0 0.033 6641.2 0.466 
Wavelength (A) 6511.1 0.053 6646.2 0.230 
Figure 9. The mean transmission profiles of the individual glass segments, 6516.1 0.086 6651.2 0.219 
A to D (cyan, blue, green and red, respectively), making up the Ha filter 6521.1 0.136 6656.3 0.133 
and the overall mean profile (black). 6526.1 0.208 6661.3 0.081 
6531.2 0.307 6666.3 0.048 
6536.2 0.429 6671.3 0.029 
while C and D span the NE and NW, respectively. Table 2 identifies 6541.2 0.575 6676.4 0.018 
the mosaic CCDs beneath each segment, and sets down the centre- 6546.3 0.700 6681.4 0.011 
to-corner range in CWL and the typical throughput integral. The 6551.3 0.799 6686.4 0.007 
laboratory tests showed us that the CWL of segments A, C and D is 6556.3 0.868 6691.5 0.005 
shortest in the segment centre, and drifts longwards according to a 6561.4 0.915 6696.5 0.003 
. . 6566.4 0.936 6701.5 0.002 
centrosymmetric pattern, as the corners and sides are approached. 
pecs p 6571.4 0.947 6706.5 0.002 
For segment B, the centre-to-corner drift is reversed, with the result 65764 0.954 67115 0.001 
that the corner CWLs are bluer than in the centre of the glass. 6581.5 0.959 6716.5 0.001 
Segment B also has the highest mean full width at half-maximum 6586.5 0.961 6721.6 0.000 


(FWHM) and highest average peak transmission: integrated over 
the bandpass this is a difference in throughput of 0.045 mag relative 
to A, C and D. The pipeline-applied illumination correction aims 
to eliminate this contrast. Area-weighted transmission profiles for 
the four segments are shown in Fig. 9, along with the overall mean 
profile. The latter is also given numerically in Table 3. 

Compared to the Ha filter used in the IPHAS survey, NB- 
659 has a CWL that is redder on average by ~20 Å, is around 
10 percent wider and has a higher overall throughput leading 
to zero-points ~0.2 higher. The known variations of bandpass 
across the four segments have implications for how best to exploit 
VPHAS+ data. To anticipate these, we have carried out two types 
of simulation based on the lab measurements in order to identify 
them. We describe these next, and summarize the implications in 
Section 3.4. 


3.2 Simulation of the main stellar locus in the (r — Ha, r — i) 
diagram 


To gain an impression of the extent of the uniformity of performance 
with regard to normal main-sequence stars, (r — Ha, r — i) tracks 
were (1) computed for each measured Ha transmission profile using 
exactly the same method as was followed by Drew et al. (2005) for 
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the analysis of IPHAS data, (2) rescaled to a common integrated 
throughput, mimicking the effect of the pipeline illumination cor- 
rection and (3) compared to the mean pattern by subtracting off the 
computed mean track. The result of this is shown in Fig. 10. The 
track differences picked out in red are from the segment corners 
exhibiting the largest CWL shifts. It can be seen that the tracks 
follow the same trend to within +0.02 up to about r — i = 1.2 
(corresponding to M3 spectral type), after which there is a clear 
fanning out. This shows that the obtained r — Ha excesses should 
fall within the target photometric precision range of the survey for 
all except mid- to late-M stars. 

The sensitivity of the M stars to variations in the narrow-band 
transmission profile is a point of note, while not actually a surprise. It 
arises from the great breadth of the feature in M-star spectra created 
by the absorbing TiO bands displaced to either side of the narrow Ha 
bandpass, and the fact that the resulting inter-band flux maximum 
falls at wavelengths shortwards of Ha. As these molecular bands 
strengthen with increasingly late spectral type, the  — Ha apparent 
excess grows along with the sensitivity to the exact placement of 
the bandpass. Viewed in these terms, the (r — Ha, r — i) colours 
of unreddened mid- to late-M dwarfs provide an empirical gauge 


€zoz 1sn6ny g, uo sanb Áq 6zP6L01/9E0Z/E/OrP/elO1UWe/SesULU/WOD' dno‘s1Wapede//:sdyYy Woy popeojumoq 


— Ha) deviation 
8 
T 


(r 
[i [i 
[e] o 
So [e] 
D A 
T T 
i i 


| 
o 
o 
& 

T 

| 


-0.10 fe fi Se ga E E E E 
-0.2 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 


r= 1 


Figure 10. The r — Ha deviations computed for all measured positions 
on the Ha filter, NB-659, after correcting the data for bandpass-integrated 
throughput variations. The data are colour-coded according to the position of 
measurement as in Fig. 8. The most discrepant corner positions, all plotted 
in red, are located in segments C and D. 


of filter bandpass uniformity and/or typical CWL. To minimize 
the bandpass sensitivity and hence spread seen at late M, the CWL 
would need to be lowered to around 6530 A or less. In practice, later 
M dwarfs are sufficiently faint that they normally appear in (r — 
Ha, r — i) diagrams as a relatively sparse distribution of scarcely 
reddened objects — falling within a thinly populated continuation of 
the unreddened main sequence, redwards of (r — i) = 1.0, rising 
from r — Ha ~ 0.5 up to ~0.8 (see Fig. 17 or 20). Reddened M 
dwarfs are usually just too faint to be detected. 

Selection of mid-to-late M dwarfs is therefore straightforward, 
but quantitative interpretation of (r — Ha) should be presumed more 
uncertain than at earlier spectral types. Similar effects will be seen 
in the M-giant spur located at lower (r — Ha) in the (r — Ha, 
r — i) diagram (see Fig. 16). However, as red giants will be picked 
up by VPHAS-+ at large distances through significant reddening, 
a precautionary check on the impact of non-zero extinction on this 
fanning in colour has been made: tracks of the type compared in 
Fig. 10 were recalculated for Ay = 6 and no noticeable additional 
effect was found (see also Fig. 17). 


3.3 Simulation of the impact of source radial velocity on 
in-band emission line fluxes 


Simulations have also been performed to consider how the filter 
captures emitted Ha flux, as a function of location within the field 
of view and source radial velocity. An ideal filter, centred on the 
mean rest wavelength of the imaged Hæ emission and placed in 
a high f-ratio optical path, would be insensitive to radial-velocity 
shifts up to a limit proportional to the FWHM of the bandpass. 
The desired capabilities of the VPHAS+ Ha filter are separation 
between Ha emission line objects and the main stellar locus — and, 
better still, a regular mapping of measured r — Ha excess on to 
emission equivalent width (cf. Drew et al. 2005, fig. 6). The two 
representative spectra used to investigate how these capabilities are 
affected by changing source radial velocity are shown in Fig. 11. 
Both spectra were blueshifted to —500 km s~!, and then shifted 
redwards in steps of 100 km s~! at a time, up to +500 km s™! (alto- 
gether a displacement of 22 A)- calculating at each step the integral 
of the spectrum folded through the filter transmission profile. The 
resultant in-band fluxes were converted to magnitudes, and then 
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Figure 11. Top panel: an example of a very bright, simple Hæ emission 
profile (taken from Corradi et al. 2010). The emission equivalent width is 
220 Å, and the FWHM of the observed profile is close to 390 km s~!. The 
mean radial velocity of the line is +35 km s~!. The difference between a 
pure continuum magnitude and that including the line is somewhat in excess 
of 1.2. Lower panel: the contrast of the line relative to continuum is much 
less here (EW ~ 20 A), and the FWHM is somewhat wider at 570 km s~! 
(a classical Be star, taken from Raddi et al. 2013). The mean radial velocity 
is —50 km s™!. The continuum-only magnitude is fainter by about 0.2 only, 
here. 


shifted by the amount required to match the integrated transmission 
to the overall mean for the filter (again mimicking the function of 
the pipeline illumination correction). In real use, we would expect 
the majority of emission line objects to present with FWHM no 
greater than either of these examples (interacting binaries and WR 
stars do present with much broader emission, however). 

The radial-velocity range explored was chosen with the follow- 
ing considerations in mind. Emission line stars in the thin disc will 
commonly have radial velocities falling within the range — 100 to 
+100 km s~!. In the bulge, larger radial velocities may be encoun- 
tered: excursions to +200 km s~! are observed in CO (Dame, Hart- 
mann & Thaddeus 2001) within ~20° of longitude of the Galactic 
Centre, and for a minority of inner-Galaxy PNe, radial velocities 
have been obtained that extend the range almost to +300 km s~! 
(see Durand, Acker & Zijlstra 1998; Beaulieu et al. 2000). 

The results of this exercise are plotted in Figs 12 and 13. 
Fig 12 shows that segments A and B come very much closer 
to independence of radial velocity than segments C and D, in 
terms of measured Ha excess. In all segments, reasonable fidelity 
(a flat or nearly flat response) is achieved around the segment 
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Figure 12. Results of simulation of the in-band flux as a function of source 
radial velocity for the EW = 220Å emission spectrum (upper panel in 
Fig. 11). The fluxes are expressed as magnitude offsets relative to the peak 
simulated in-band flux. The panels representing the segments are arranged 
as they are imposed on the plane of the sky, i.e. A covers the SW quadrant 
and D covers the NW when the filter is stored in OmegaCAM’s magazine 
B. The curves are colour-coded according to measurement position as in 
Fig. 8. Galactic sources will usually fall well within the marked velocity 
range, —200 to +200 km s~!. The more problematic red curves, represent- 
ing the response of the segment corners, account for ~8 percent of each 
segment’s unvignetted area only. 


0.0 + + 


0.1 F > 


0.27 + 


0.0 + + 4 


Ha magnitude offset 


| 


0.1 b + 4 


0.2 + + | 


1 1 L E 1 L l 1 L js 1 l 1 L 1 1 L 1 1 L L 1 
—400 -200 0 200 400 —400 -200 0 200 400 
velocity offset (km/s) 


Figure 13. Results for an emission line net equivalent width of 
20 Å(spectrum shown in the lower panel of Fig. 11). Otherwise as Fig. 12. 


centre — although in B, uniquely, the corners happen to perform 
a little better than the centre. Clearly segment D, where the trans- 
mission is centred on longer wavelengths than in the other segments, 
could yield measurements in its corners (perhaps 8 per cent of its 
unvignetted area) of Hæ magnitude, or of r — Hø, that underesti- 
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Figure 14. The layout of the Hæ filter, showing the positions of the PN- 
test measurements, colour-coded according to flux relative to the corners of 
segment B. Those shown as black crossed boxes have a scaled flux >0.95. 
Those in blue are a little lower with a scaled flux in the range 0.90-0.94. 
The two positions marked in green have fluxes scaling to 0.88 and 0.89, 
while the two in red scale to 0.84 and 0.77 (both in segment D). There is no 
measurement for the centre of segment C because of a telescope pointing 
imprecision. The other centre measurements were obtained ~2 arcmin away 
from the true centres plotted in order to avoid CCD gaps. The dashed lines 
mark the limits of the T-bar vignetting. 


mate the flux of similarly high-equivalent-width Ha emission by 
up to ~0.3 (out of a true excess, expressed in magnitudes of ~1.2). 
Segment C performs similarly, but the potential flux drop associated 
with its corners is less pronounced. 

We can compare the expectations created by these simulations 
with the results of an on-sky experiment in which the PN ESO 
178-5 (or PNG 327.1—02.2) has been exposed in Ha at a series 
of positions in the image plane, placing it well into every corner 
and also close to the centre of each of the four filter segments 
(see Fig. 14). The observed integrated counts variation might be 
predicted to be somewhat stronger than in Fig. 12 given that the in- 
band continuum flux from this PN will be relatively even weaker. 
But this will be offset by the additional flux due, in particular, to 
[N u] 46584. The PN chosen for this test was picked both because 
it is well calibrated (Dopita & Hua 1997) and because its local 
standard of rest (LSR) radial velocity is quite large and negative 
(U_sp = —88.7 km s~!: on 2013 April 19 when it was observed, 
this will have shifted to —105 km s~! at the telescope). It also 
happens to possess [N 11] 46584 emission that is scarcely less bright 
than Ha (the former has 98 percent the flux of the latter: Dopita 
and Hua also provide a spectrum of this nebula). 

Background-subtracted aperture photometry of the PN and a 
moderately bright star nearby, serving as a continuum reference, 
was carried out on the reduced images. These measurements reveal 
a pattern of behaviour that essentially tracks the results shown in 
Fig. 12: the continuum reference itself shows a total count variation 
of +5 percent across all pointings, while the PN counts, after scal- 
ing to the reference, range from +5 percent down to —23 per cent 
relative to the values for the corners of segment B. In the extreme 
case that all of the Ha and [N 1] 46548 emission had been shifted out 
of the bandpass, the maximum drop for this PN would be 
—57 per cent (the remaining 43 per cent being attributable to [N 1] 
46584). The pattern across the filter of the results emerging from 
this trial is shown in Fig. 14. 
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The radial-velocity dependence of the transmitted flux may ac- 
cordingly become an issue for objects with very strong Ha emission 
where the aim is accurate flux determination, unless attention is paid 
to where the object falls in the image plane. Qualitatively, the issue 
is less critical: regardless of where the object is located, the changes 
in transmission are not so large that there will be frequent failures 
to distinguish strong Ha emitters — i.e. they will still appear above 
the main stellar locus in the (r — Ha, r — i) diagram. In the ex- 
ample simulated, the outer reaches of segment D would bring (r — 
Ha) down to 0.9-1.0, a level that nevertheless remains clear of the 
domain that might be occupied by unreddened, non-emission-line 
very late type M dwarfs [(r — Ha) ~ 0.8), cf. Fig. 16]. As further 
context, we note that nearly continuum-free emission line objects, 
such as PNe and H 1 regions, present with r — Ha ~ 3. 

Where the line emission itself contributes only a minority of the 
measured narrow-band Ha flux, the trends seen are much more 
subdued (Fig. 13). For the example shown, only 20 per cent of the 
total in-band flux is attributable to the net line emission, rather than 
most of it as in Fig. 12. Again, the corners of segments C and 
D perform least well, in under-representing the emission flux by 
up to ~0.05 mag at the most negative likely Galactic plane radial 
velocities. Otherwise, the performance is predicted to be within the 
anticipated 0.02-0.03 error budget of the survey. 


3.4 Implications of the Ha filter properties for VPHAS+ and 
its exploitation 


In summary, for most purposes the Ha filter performs as required, 
and has very good throughput. For the great majority of stars making 
up the main stellar locus, there will be the desired fidelity of the 
(r — Ha) colour, and the great majority of emission line objects will 
be detected with the same facility as they are by IPHAS. 

There are two caveats to note. First, in Section 3.2 it was shown 
that variations in CWL across the filter segments will lead to thicken- 
ing of the loci traced by mid-to-late M stars. These same variations, 
of what is a relatively red Ha passband, also introduce the potential 
for underdetermination of Ha fluxes for objects/nebulosity in parts 
of the image plane for sources with significantly negative radial 
velocities (Section 3.3). This becomes most serious for emission 
line sources falling near the vignetted corners of segments D and 
C, where a 20—30 percent undercounting in Ha may occur for ra- 
dial velocities approaching —200 km s™!. As is always the case for 
narrow-band Ha filters, the common presence of significant [N 11] 
226548, 6584 emission bracketing Ha in PNe or H 1 regions com- 
plicates the expected signature. However, it can be guaranteed in all 
but rare, exotic circumstances that the stronger 46584 component 
of the [N u] doublet falls well within the bandpass. 

The outstanding practical consequence of the filter’s transmis- 
sion characteristics for the survey strategy is that, for quantitative 
reliability, measures of (r — Ha) obtained using segments A and 
B, and the central zones of C and D (out to ~13 arcmin) are to be 
favoured. This appreciation is half of the reason for the adoption of 
offsets of several arcminutes between the three successive pointings 
made in this filter (Fig. 3) for each field — our strategy ensures that 
objects captured in segment corners in the first pointing fall close 
to segment centres in the third. 

The rest of the motive is to mitigate the cross-shaped vignetting 
due to the blackened T-bars holding the segments in place (Fig. 7). 
Each arm of the cross casts a shadow entirely contained within a 
strip 4 arcmin wide. By choosing to offset at least this much in both 
RA and Dec. between the three exposures obtained per field, we 
raise the probability of at least one high-confidence r — Ha colour 
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measurement per detected source in the final catalogue to very 
nearly 100 percent, and the probability of two to over 95 per cent. 

Finally, we remark that the combination of large offsets and three 
pointings has the consequence that the fraction of sky within the 
survey footprint missed altogether, due to dead areas between the 
CCDs and vignetting, is under 0.3 per cent. 


4 SIMULATION OF VPHAS+ STELLAR 
COLOURS 


The five bandpasses of the survey provide the basis for the construc- 
tion of a range of magnitude—colour and colour—colour diagrams. To 
take full advantage of them, knowledge is needed of the behaviours 
that can be expected of the colours of normal stars. 

We have simulated colours for solar-metallicity main-sequence 
and giant stars using the same method as employed by Sale et al. 
(2009). We adopt the definition of these two sequences in (Teff, 
log g) space given by Straizys & Kuriliene (1981). Then for each 
spectral type along these sequences, solar-metallicity model spectra 
were drawn from the Munari et al. (2005) library. At a binning of 
1 A, the spectra in this library are well enough sampled to permit the 
calculation of narrow-band Ha relative magnitudes with confidence, 
alongside the analogous broad-band quantities. More detail on the 
broad-band filter transmission profiles, shown in Fig. 2, and on the 
CCD response curve is provided on the ESO website.* To ensure 
compliance with the Vega-based zero-magnitude scale, we have 
defined the synthetic colour arising from a flux distribution F, as 
follows: 


[TAF dA 


[TAFA 
mı — m = —2.5 log or ; 
[TAF van 


2.5 log | 2 
|+ 08 | tr dk 


where Tı and T are the numerical transmission profiles for fil- 
ters 1 and 2, after multiplying them through by the atmospheric 
transmission (Patat et al. 2011) and mean OmegaCAM CCD re- 
sponse curves. The SED adopted for Vega, Fy.y, is that due to 
Kurucz (http://kurucz.harvard.edu/stars/VEGA/), Tere = 9550 K, 
log g = 3.95. Where needed for comparison, we have also com- 
puted colours based on the Pickles (1998, hereafter P98) spec- 
trophotometric stellar library [the approach adopted by Drew et al. 
(2005) for IPHAS]. To maintain precision, the numerical quadrature 
resamples the more smoothly varying transmission data on to the 
sampling interval of the stellar SED. 

The (r — Ha) excess is evaluated in exactly the same way as 
the broad-band colours. Since Vega is an AOV star, its SED at Ha 
incorporates a strong absorption line feature that reduces the in- 
band flux below the pure continuum value. Unlike the broad-bands, 
the Ha narrow band has not yet been standardized and so there is not 
a formally recognized flux scale. However, we can specify here that 
the integrated in-band energy flux for Vega, on adopting the mean 
profile for the VST filter, is 1.84 x 1077 ergs cm~? s7! (at the top of 
the Earth’s atmosphere). To assure zero colour relative to the optical 
broad-bands, this flux is required to correspond to mya ~ 0.03. The 
reduction in zero-point (zpt) that the computed in-band flux implies 
relative to the flux captured by the much broader r band — based on 
folding Vega’s SED with lab measurements of the filter throughputs 
corrected for atmosphere and detector quantum efficiency — is 3.01. 
Current practice in VPHAS-+ photometric calibration is accordingly 
to adopt zpt(NB-659) = zpt(r) — 3.01 mag as the default calibration 
for the narrow band: in Section 5 where a direct comparison is 


3 http://www.eso.org/sci/facilities/paranal/instruments/omegacam/doc. 
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made with SDSS spectroscopy, this offset is found to be satisfactory. 
When applied, it assures that data obtained in photometric, or stable, 
conditions yield zero r — Ha colour for AO stars. 

Both main-sequence and giant-star colours have been calculated 
for a range of reddenings and optical—IR extinction laws as for- 
mulated by Fitzpatrick & Massa (2007). The unreddened colours 
for the mean Galactic law (R = 3.1) are laid out here in Table 4. 
Appendix A provides additional tables that specify the colours of 
main-sequence stars at selected reddenings and for two further rep- 
resentative reddening laws (R = 2.5 and 3.8). These can be used 
to construct intrinsic-colour-specific reddening lines. For the large 
range in extinction sampled along many Galactic plane sightlines, 
these reddening trends are slightly curved (see the examples shown 
ine.g. Sale et al. 2009). In this paper, we use Ao, the monochromatic 
reddening at 5500 A to parametrize the amount of reddening, rather 
than the band-averaged measure, Ay. In most circumstances, these 
quantities are almost identical. 


Based on the data from these tables, the main-sequence and giant 
tracks are as shown in Figs 15 and 16. These identify where the main 
stellar loci will fall. It is important to note that the OmegaCAM u 
filter, like all filters constructed for this challenging band, exhibits 
a low-level red leak. In this instance, lab measurements indicate 
transmission at levels between 1075 and 10~* within limited win- 
dows around ~9000 A. This is enough to begin to noticeably, and 
erroneously, brighten the u magnitudes of normal stars reddened to 
g —r> 3. Because of this, and because the measurement of very low 
level leakage is itself subject to proportionately higher uncertainty, 
we do not plot or tabulate u — g data beyond the g — r = 3 limit. 
Very few detected sources are so extreme. In practice, VPHAS+ 
u — g is faithful for extinctions up to Ap ~ 6, but gradually thereafter 
it transforms into a colour that behaves crudely as —(g — z). 

The (r — Ha, r — i) colour—colour diagram is not subject to such 
effects, and therefore remains sound across a wider spread in visual 
extinctions. Synthetic tracks are presented in Fig. 16 for Ag = 0, 2, 


Table 4. VST/OmegaCAM synthetic colours for unreddened main-sequence dwarfs and giants. 


Spectral type Main sequence (V) Giants (IIT) 
(model) (P98 library spectra) 

u— g g-r r—i r — Ha u— g g-r r-i r — Ha r-i r— Ha 
06 1.494 0.313 0.145 0.071 
08 1.463 0.299 0.152 0.055 —0.158 0.074 
09 1.426 0.271 0.142 0.064 1.426 0.271 0.142 0.064 
BO 1.404 0.267 0.143 0.058 1.404 0.267 0.143 0.058 
Bl 1.296 0.236 0.130 0.052 1.316 0.234 0.130 0.057 —0.095 0.071 
B2 1.181 0.214 0.117 0.049 1.209 0.211 0.116 0.056 
B3 1.025 0.182 0.098 0.048 1.046 0.182 0.098 0.054 —0.035 0.083 
BS 0.799 0.133 0.071 0.043 0.814 0.134 0.072 0.050 —0.016 0.083 
B6 0.699 0.116 0.062 0.040 0.714 0.116 0.062 0.046 
B7 0.550 0.094 0.051 0.033 0.568 0.095 0.051 0.041 
B8 0.361 0.071 0.039 0.022 0.383 0.072 0.039 0.032 
B9 0.168 0.040 0.023 0.009 0.186 0.044 0.024 0.021 —0.018 0.035 
AO —0.024 0.000 —0.003 —0.002 0.030 0.009 0.006 0.011 0.012 0.034 
Al 0.007 0.015 0.004 —0.004 0.007 0.004 0.000 0.008 
A2 0.039 0.038 0.014 —0.005 0.051 0.022 0.009 0.008 
A3 0.064 0.062 0.025 —0.005 0.085 0.045 0.019 0.007 0.037 0.063 
AS 0.096 0.130 0.056 0.008 0.143 0.107 0.048 0.015 0.096 0.087 
A7 0.073 0.206 0.089 0.030 0.145 0.179 0.078 0.032 0.115 0.094 
FO 0.003 0.336 0.153 0.086 0.091 0.317 0.144 0.084 0.156 0.102 
F2 —0.021 0.396 0.182 0.111 0.064 0.380 0.174 0.109 0.204 0.172 
F5 —0.039 0.505 0.230 0.150 0.046 0.491 0.224 0.148 0.238 0.164 
F8 —0.013 0.587 0.263 0.174 
G0 0.012 0.628 0.278 0.185 0.333 0.213 
G2 0.011 0.628 0.279 0.185 0.253 0.759 0.329 0.215 
G5 0.162 0.756 0.327 0.217 0.405 0.870 0.368 0.235 0.396 0.250 
G8 0.355 0.845 0.358 0.233 0.531 0.944 0.395 0.247 0.419 0.247 
KO 0.523 0.938 0.396 0.248 0.640 1.002 0.417 0.256 0.446 0.254 
Kl 0.551 0.954 0.403 0.251 0.803 1.085 0.451 0.269 0.468 0.269 
K2 0.629 0.993 0.419 0.258 0.963 1.159 0.483 0.281 0.508 0.293 
K3 0.779 1.062 0.447 0.269 1.227 1.276 0.536 0.299 0.514 0.286 
K4 0.871 1.108 0.468 0.278 1.374 1.342 0.585 0.320 0.592 0.313 
K5 1.083 1.210 0.522 0.300 1.578 1.420 0.630 0.338 0.714 0.337 
K7 1.387 1.402 0.724 0.387 
MO 1.372 1.411 0.789 0.411 1.697 1.454 0.686 0.360 0.827 0.411 
MI 1.335 1.439 0.934 0.467 1.838 1.506 0.769 0.395 0.872 0.401 
M2 1.262 1.442 1.112 0.522 1.938 1.546 0.825 0.420 0.920 0.443 
M3 1.236 1.447 1.179 0.545 1.980 1:536 0.933 0.444 1.165 0.471 
M4 1.248 1.457 1.168 0.543 1.959 1.551 1.136 0.500 1.472 0.512 
M5 2.009 1.569 1.296 0.556 1.739 0.560 
M6 2.199 1.612 1.267 0.554 
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Figure 15. The expected positions of main-sequence and giant stars in the 
(u — g, g — r) plane. For the main sequence, tracks are shown for the 
monochromatic reddenings Ag = 0, 2, 4, 6 and 8 (working from left to 
right). The red leak in the u filter starts to lower u — g, noticeably from 
Ao = 6 (tracks drawn in red). The giant-star tracks, drawn as dashed lines 
for Ap = 0, 2 and 4 only, are very similar to their main-sequence counterparts 
except at the latest types. 


r — Ha 


Figure 16. The expected position of main-sequence and giant stars in the 
(r — Ha, r — i) plane. Tracks are shown for the monochromatic reddenings 
Ao = 0, 2, 4, 6, 8 and 10 (from left to right). Solid lines represent the main- 
sequence tracks, while dashed lines are used for the giant tracks. The lines 
in red are giant-star tracks derived from P98 spectrophotometry. 


4, 6, 8 and 10. The main-sequence tracks shown are similar to those 
appropriate to IPHAS (cf. fig. 6 of Sale et al. 2009). But a problem 
emerges when it comes to the simulation of red-giant colours. Purely 
theoretical simulation predicts late-K and M giant colours closely 
resembling those of dwarfs, whereas simulation using P98 library 
spectrophotometry indicates a distinctive flattening of the M-giant 
track, peeling away from the steadily rising main-sequence track. 
Fig. 16 points out this contrast. Inspection of Table 4 reveals this 
is a problem linked mainly to simulation of the i spectral range, 
which renders (r — i) progressively larger from late-K into the M 
giant range when the library spectra are used in place of model 
atmospheres. 

Evidence that M giants are better reproduced by synthetic pho- 
tometry based on flux-calibrated spectra is provided by Fig. 17. This 
figure also compares new VPHAS-+ data with their cross-matches 
in the IPHAS survey within an ~0.2-square-degree equatorial field 
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Figure 17. Equatorial VPHAS-+ data (upper panel) and IPHAS data (lower 
panel) compared to show the different appearance of M-giant (r — Ha, r — 
i) colours. The photometry is extracted from an ~0.2-square-degree region 
of sky, centred on £ = 35°95, b = —3°13. The magnitude range is limited 
to 13 < r < 18 in both cases. Telescope-appropriate giant tracks, computed 
from the P98 spectrophotometric library, for Ag = 2 and 4 are superimposed 
in red in both panels. In the upper panel, giant tracks computed from model 
atmospheres for the same reddenings are shown in blue. The black dashed 
line in each panel shows the synthetic early-A reddening line, from Ag = 0 
to 10. The stars appearing in the IPHAS catalogue as candidate emission 
line stars are enclosed in cyan boxes in both panels. 


(£ = 35°95, b = —3°13), and shows selected synthetic tracks super- 
imposed. The photometry from the two surveys of the most densely 
populated part of the main stellar locus to (r — i) ~ 1.5 substantially 
overlap, but not perfectly — the response functions describing the 
three bandpasses involved in Fig. 17 undoubtedly differ in detail 
between the two telescopes. On cross-correlating either (r — i) or 
(r — Ha) between the two surveys, it becomes clear that the IPHAS 
colour has the somewhat larger dynamic range. This is the reason 
for the slightly more stretched appearance of both the main locus 
and the early-A reddening line in the IPHAS diagram relative to 
that for VPHAS+. 

At (r — i) > 1.5 in Fig. 17, it can be seen that the M-giant spurs 
look very different. First, the VPHAS+ M giants fall into a nearly 
flat distribution lying at lower (r — Ha), compared to the more 
steeply rising higher IPHAS M-giant sequence. However, as long 
as the data are interpreted with reference to telescope-appropriate 
synthetic photometry, the two data sets will lead to the same in- 
ference. In the example shown in Fig. 17, the comparisons with 
suitable synthetic giant tracks indicate that the maximal extinction 
in the field can be no more than Ag ~ 4. The extinction mea- 
sures due to Marshall et al. (2006), based on 2MASS red-giant 
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photometry, indicate a maximum Galactic extinction of Ax ~ 0.3 
for this pointing. For a typical Galactic R = 3.1 reddening law, this 
scales up to Ap ~ 3.3 (roughly — see Fitzpatrick & Massa 2009). 
If model-atmosphere giant tracks are referred to instead, the M gi- 
ants would have to be read as demanding visual extinctions ranging 
from ~4 upwards. 

Fig. 17 also demonstrates the broadening in (r — Ha) of the 
VPHAS+ M-giant sequence that was foretold in Section 3. The 
IPHAS counterpart is evidently much sharper, as it rises to higher 
r — Ha with increasing r — i. The main practical impact of this 
difference is that IPHAS M-giant photometry is the better starting 
point for picking apart chemistry differences (Wright et al. 2008). 
But it is as true of the VPHAS-+ (r — Ha, r — i) diagram as it is of 
its IPHAS equivalent — that M giants at r — i > 1.5 sit below and 
apart from M dwarfs. 

Fortuitously, Fig. 17 identifies an advantage of the generally good 
seeing available at the VST. There are two candidate emission line 
objects apparent in the IPHAS selection (enclosed in cyan boxes, in 
Fig. 17), that drop back into the main stellar locus in the VPHAS+ 
data. Inspection of the images shows that both stars are in close 
doubles of similar brightness, of under 2 arcsec separation. Because 
they are a little better resolved in VPHAS+ (~0.8 arcsec seeing), 
than in IPHAS (~0.9 arcsec seeing), the pipeline makes a better 
job of the assigning magnitudes in the different bands to the blend 
components. This example nicely illustrates the most common rea- 
son for bogus candidate emission line stars in either VPHAS+ or 
IPHAS - improperly disentangled blends. Candidate emission line 
stars should always be checked for this kind of problem before a 
spectroscopic follow-up. Otherwise, experience with IPHAS gives 
confidence that the selection of emission line objects via VPHAS+ 
will be highly efficient (see e.g. Vink et al. 2008; Raddi et al. 2013). 

Finally, it is worth noting that the bright limit of the survey at 12— 
13th magnitude effectively excludes any unreddened stars of earlier 
spectral type than ~GO. Before more luminous stars of spectral 
type F and earlier can enter the survey sensitivity range, they need 
to be at distances in excess of 1 kpc, typically, where low extinction 
becomes increasingly improbable. This constraint bestows a signif- 
icant selection benefit in that only unreddened or lightly reddened 
subluminous objects, with intrinsically blue colours, are left stand- 
ing clear near the blue end of the main stellar locus in commonly 
constructed photometric diagrams. In this domain, VPHAS+ has 
important selection work to do. 


5 PHOTOMETRY VALIDATION: A 
COMPARISON OF SDSS AND VST DATA 


Before the start of survey field acquisition, we obtained observations 
in all survey filters of two pointings that fall within the SDSS pho- 
tometric and spectroscopic coverage (Abazajian et al. 2009). These 
were centred on RA 20°47™ :5387, Dec. —06°04' 145 (J2000) and 
RA 21° 04™ 25894, Dec. +00° 59’ 15”8 (J2000) — fields that happen 
to include a number of white dwarfs and cataclysmic variables (not 
discussed further here). The main aim of the data was to verify 
VST photometry both by comparison to SDSS photometry and to 
synthetic photometry derived from SDSS spectra. The VST obser- 
vations were obtained on 2011 September 21, during clear weather 
at a time of generally sub-arcsecond seeing. The exposure times 
differ only a little from those now in general survey use: the g ex- 
posures were 30 s, rather than 40, and į was exposed for 20 s rather 
than 25 — the other times were as given in Section 2.1. 

The photometry on the sources in these fields have been pipeline 
extracted and calibrated in the standard way, and have been cross- 
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Figure 18. The distributions, by band pass, of the measured magnitude 
differences between VST test data and SDSS, after correction of the latter 
magnitudes from the AB to the Vega zero-magnitude scale. These were 
obtained from two VST pointings, away from the Galactic plane, that overlap 
the SDSS footprint. The biggest differences are found in the u band where 
the VST data are fainter than the corrected SDSS values by ~0.12, in the 
median. 


matched to their SDSS counterparts. The number of cross-matched 
stars used in this exercise ranges from ~2500 (u) up to ~10000 
(r). For a star to be included, it must be unvignetted, should have 
a star-like point spread function, should lie within 0.5 arcsec of 
its SDSS counterpart and should fall within the magnitude range 
16 > rysr > 19. The SDSS selection constraints were set to exclude 
blended and saturated sources, and sources close to detector edges. 
In addition, it was required of every source that, in both surveys, 
the formal error on the magnitude measurement is less than 0.03. 

In Fig. 18, we plot the histograms of the magnitude differences 
between the two surveys, according to passband — pooling the data 
from both pointings. If the starting assumption is that the VST 
broad-band filters are identical to the SDSS set, the predicted mag- 
nitude for each star in each filter is the measured SDSS magnitude, 
less the offset between the AB and Vega scales (essentially the nu- 
merical difference between the magnitude of Vega in the AB system 
and its value of 0.02—0.03, according to the alternative Vega-based 
convention — see table 8 in Fukugita et al. 1996). In the g, r and 
i bands, the predicted and observed magnitudes are well enough 
aligned, and the interquartile spread is consistent with the way the 
data were selected for random errors less than 0.03. However, in u 
there is a discrepancy that exceeds the expected error: the median 
difference deviates by 0.12 mag and the width of the distribution is 
twice that arising in the comparisons of the other bands. 

The fuller picture is presented in Fig. 19 which shows the broad- 
band magnitude differences as a function of the relevant SDSS 
colour for the second of our two fields (only). In all four passbands, 
including u, the colour dependence can be seen to be very weak 
in that the loci traced out by the plotted stars are — to a first ap- 
proximation — flat. The discrepancy seen in u is revealed as mainly 
a zero-point shift, combined with scatter that exceeds the formal 
errors. In highly reddened Galactic plane fields, the stellar colour 
effects may become more pronounced as extinction modifies the 
effective sampling of the passbands. 
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Figure 19. Measured magnitude shifts between VST and SDSS cross- 
matched objects as a function of the SDSS colour, for the second of the two 
fields observed (only). The data in blue are Au versus (u — g)spss; the data 
in green and red are, respectively, Ag and Ar versus (g — r)spss; the data 
in black are Ai versus (r — i)spss. The horizontal lines show where the four 
loci would be expected to lie in the case that the SDSS and VST broad-band 
filters were identical, and the calibrations perfect. The data fall into loci that 
are not far displaced from these horizontal lines and are almost as flat: this 
indicates that the colour-dependent terms that would be needed in equations 
to transform between the SDSS and VST systems are small. 


Table 5. Mean magnitude offsets between VST pho- 
tometry and cross-matching synthetic photometry de- 
rived from the SDSS data base of spectra. The synthetic 
magnitude scale adopts magnitudes for Vega itself of 
0.026 (see Bohlin & Gilliland 2004). 


Offset Field 1 (117 stars) Field 2 (50 stars) 
uysT — Usyn 0.07 + 0.39 0.11 £0.17 
8VST — Ssyn 0.06 + 0.09 0.06 + 0.04 
TYST — "syn 0.01 0.06 0.03 0.03 
ivst — tsyn —0.06 + 0.08 —0.05 + 0.03 


AS a Separate exercise, we have used SDSS spectra to synthesize 
magnitudes and colours for stars with cross-matching VST pho- 
tometry. The spectral type range present within this much smaller 
sample runs from B-type through to early M-type (M1). At wave- 
lengths below 3800 A falling within the u band, it was necessary to 
extrapolate the spectra using appropriately chosen P98 library data. 
The result of this comparison is agreement between the VST and 
synthesized magnitudes at the ~5 per cent level (Table 5), with the u 
band as the outlier exhibiting much more pronounced scatter as well 
as somewhat higher offset. This pattern echoes the behaviour ap- 
parent in the VST-SDSS purely photometric comparison of Fig. 18, 
using a much larger sample. The difficulty is not confined to VST 
u however, in that SDSS u photometry fares scarcely any better 
relative to synthesis from the spectra (for the two fields, offset and 
scatter are —0.09 + 0.37 and —0.04 + 0.16). As more blue survey 
data are accumulated, it may become clear that the u zero-point will 
benefit from being tied to that of g for those fields observed in the 
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best conditions, as is presently done for Ha with respect to r. This 
option is not yet enacted. For the time being, it must be acknowl- 
edged that pipeline u calibration is more approximate than those of 
the other bands. 

We have also used the reduced cross-match sample to look at how 
the VST photometric r — Ha colour compares with its counterpart 
synthesized from spectra — looking, in particular, for any trends as 
a function of the distance from the field centre. No such trend is ap- 
parent, thereby meeting the expectation that the narrow-band fluxes 
of normal stars, to early-M spectral type, would not be affected 
by the pattern of bandpass shifts discussed earlier in Section 3 (cf. 
Fig. 10). However, we do find that in order to make this detailed 
comparison, systematic offsets had to be removed from the VST 
photometry first. These were 0.075 in (r — i), in the sense that the 
VST colours were too red by this amount, and a 0.02 reduction 
in (r — Ha). The (r — i) offset is consistent with the broad-band 
magnitude offsets listed in Table 5 and hence is as expected. The 
(r — Ha) adjustment is small enough (i.e. within the fit error) that 
it supports the zero-point shift of 3.01 mag between the r and Ha 
bands that was identified in Section 4. Once these colour offsets 
are applied to the VST data, the rms scatter of the photometric 
(r — Ha) colour relative to its synthetic counterpart is 0.04 for 
objects brighter than r = 19. 


6 AN EXAMPLE OF POINT-SOURCE 
PHOTOMETRY DERIVED FROM A VPHAS+ 
FIELD 


The extracted point-source photometry from the VST square-degree 
field is one of the two main data products from VPHAS-+ — the other 
being the images themselves, considered below in Section 7. We 
present an example of the two essential colour—colour diagrams in 
Fig. 20, in which band-merged stellar photometry for field 1679 is 
compared with the primary diagnostic synthetic tracks presented in 
Section 4. This field includes the sky area from which data were 
taken to construct Fig. 6 illustrating typical errors. The massive 
open cluster, Westerlund 2, is located in the NE of these pointings, 
and the field as a whole includes moderate levels of diffuse, com- 
plex Hu emission. The data presented are drawn from a sky area 
centred on RA 10" 24™ 49°, Dec. —57° 58’ 00” (J2000) that spans 
1.3 square degrees — the total footprint occupied by the two offset 
positions. 

Only objects with stellar point-spread functions in g, r and i, 
brighter than g = 20, are included in Fig. 20. Where two sets of 
magnitudes are available, the mean values have been computed and 
used. A further requirement imposed is that the random error in all 
bands may not exceed 0.1. The same ~37 000 objects are included 
in both diagrams. In order to obtain the diagrams shown, the pipeline 
photometric calibration was checked and refined as follows: we 


(i) cross-matched brighter stars to APASS g, r and i photometry; 

(ii) computed the median magnitude offset (applying no colour 
corrections — it was shown in Fig. 19 that these are modest); 

Gii) corrected all g, r and i for these offsets; 

(iv) corrected u by determining the vertical shift needed in the 
(u — g, g — r) diagram to align the main stellar locus with the 
unreddened main sequence and the GOV reddening line; 

(v) corrected the Hæ zero-point and hence all Ha magnitudes 
according to the requirement that zpt(Ha) = zpt(r) — 3.01. 


This resulted in the following broad-band corrections: 
Ai = —0.004, Ar = —0.032 (red-filter set), Ar = —0.033 
(blue-filter set), Ag = 0.069 and Au = —0.31. As expected, the 
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u—g 


Figure 20. The left-hand panel shows the (u — g, g — r) diagram, the right the (r — Ha, r — i) diagram pertaining to VPHAS+ survey field 1679. Both 
diagrams are plotted as two-dimensional stellar-density histograms, rainbow colour-coded such that high source densities (80-90 per bin) are red and the lowest 
densities (one per bin) are dark blue. The binning is 0.017 x0.025 in the left-hand panel and 0.013 x 0.008 in the right. The synthetic unreddened main sequence 
is drawn in, in black, in both panels. The GOV and A3V reddening lines obtained for R = 3.1, drawn as black dashed lines, are included in, respectively, the 


(u — g, g — r) and (r — Ha, r — i) diagrams as useful aids to interpretation. 


correction that had to be applied to the u photometry was, by far, 
the largest. 

The main stellar locus can be seen to be tightly concentrated in 
both the blue and the red diagrams, and to favour lightly reddened 
G and K stars. The superimposed synthetic reddening lines [GOV 
in the (u — g, g — r) diagram and A3V in (r — Ha, r — i)] have been 
drawn adopting the R = 3.1 reddening law widely regarded as the 
Galactic norm. The blue diagram provides examples of three distinct 
typical populations falling outside the main stellar locus. Below it, 
at (u — g) > 1.5 and (g — r) > 1.5 (roughly) the plotted objects 
will mainly be M giants. Above the main stellar locus towards the 
red end, in the ranges 0 < (u — g) < 0.5 and 1.5 < (g — r) < 2.0 
lie the OB stars in and around Westerlund 2. Finally, the modest 
scatter of blue objects lying above the GOV line roughly in the 0 < 
(g — r) < 1 range will include intrinsically blue lightly reddened 
subluminous objects. 

It is interesting to note in the red diagram that there is some 
evidence that early-A stars making up the lower edge of the main 
stellar locus would better follow a different law, with R ~ 3.8 (see 
the tables in Appendix A). Indeed, a reddening law of this type has 
been inferred for the OB stars in Westerlund 2 by Vargas Alvarez 
et al. (2013). Most of the thin scatter of points below the main 
stellar locus, and some of the scatter above, in this same diagram 
will be the product of inaccurate background subtraction in Ha. 
But many of the objects lying above the main stellar locus will 
indeed be emission line objects, and some of the stars below will be 
white dwarfs. As expected, the red spurs of M dwarfs and M giants 
are broader features than their IPHAS counterparts (cf. Fig. 17 and 
associated remarks). 

For more discussion of these colour—colour diagrams, the reader 
is referred to Groot et al. (2009, UVEX) and Drew et al. (2005, 
IPHAS). 
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7 NEBULAR ASTROPHYSICS WITH VPHAS+ 
IMAGES 


Just over a decade ago, the SuperCOSMOS Ha Survey (SHS; Parker 
et al. 2005) had only just completed. This was the last survey using 
photographic emulsions that the UK Schmidt Telescope undertook. 
The 3 h narrow-band Ha filter exposures reach a very similar lim- 
iting surface brightness to the 2 min exposures VPHAS-+ is built 
around. Hence, the differences in capability are not about sensi- 
tivity, as this is roughly the same in the two surveys. Instead, it 
is about the great improvement in dynamic range on switching to 
digital detectors, the good seeing of the VST’s Paranal site and the 
added broad-bands. 

SHS, with its enormous 5° diameter field, has been comprehen- 
sively trawled for southern PNe (the Macquarie/AAO/Strasbourg 
Ha or MASH catalogue; Parker et al. 2006; Miszalski et al. 2008). 
The remaining discovery space for resolved nebulae is expected to 
be at low surface brightnesses in locations of high stellar density, 
and in the compact domain around and below the limits of the typical 
spatial resolution of SHS (~0.5-3.0 arcsec). Both these conditions 
will most often be met in the Galactic bulge, at a mean distance 
of ~8 kpc. Data taking in the bulge and its maximally dense star 
fields is planned to begin in mid-2014. 

Among PNe, small angular size is due either to great distance or 
to youth — the study of either compact category provides exciting 
possibilities. As well as the bulge, the less studied outer parts of the 
Galactic plane should be searched. In this respect, IPHAS, with its 
direct view to the Galactic anticentre, is better positioned: the on- 
going study of the anticentre PN population has revealed dozens of 
new candidates (Viironen et al. 2009a), including the PN with the 
largest galactocentric distance to date (20.8 + 3.8 kpc; Viironen 
et al. 2011). By following up such finds to measure chemical 
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Figure 21. Two PNe, NGC 2438 (top) and NGC 2899 (bottom), as they appear in the SHS and VPHAS-+ surveys. The SHS images are shown in the left-hand 
panels, with the VPHAS-+ images to the right. The bands used to form them are: NGC 2438 — SHS R/G/B = Ha/SR/SSS Bj, VPHAS+ R/G/B = Ha/r/i, NGC 
2899 — SHS R/G/B = Ha/SR/SSS Bj, VPHAS+ R/G/B = Ha/r/g. The cut-out image dimensions are 300x300 arcsec” for NGC 2438 and 200x 180 arcsec” 


for NGC 2899. 


abundances, crucial beacons are obtained for the study of the Galac- 
tic abundance gradient and its much disputed flattening towards the 
largest galactocentric radii. VPHAS+ completed the access to the 
outer plane over the longitude range 215° < £ < 270°. 

Data from both IPHAS and VPHAS+ can make fundamental 
contributions to the study of very young PNe — particularly by help- 
ing to solve the two-decades-old puzzle of how PNe already emerge 
with the observed wide variety of morphologies (round, elliptical, 
bipolar, multipolar, point-symmetric, etc. — see Sahai, Morris & 
Villar 2011). What does this variety say about the properties of 
their asymptotic giant branch (AGB) progenitors? Detailed studies 
of objects in the phases preceding the PN phase — AGB and post- 
AGB stars, proto-PNe and transition or PN-nascent objects — are 
underway (e.g. Sanchez-Contreras & Sahai 2012). Superb imag- 
ing capabilities like those of the VST, accessed via VPHAS-+, will 
support this work. 

Indeed, there is a serious paucity of very small PNe in the ex- 
isting optical catalogues: there are no PNe with angular extent less 
than 3 arcsec in the MASH catalogue (out of 903 objects; Parker 
et al. 2006), and only 8 PNe in the catalogue by Tylenda et al. 
(2003, 312 objects) in the size range 1.4—3 arcsec. There is just 
one with a confidently measured diameter below 1 arcsec in the 


larger Strasbourg Catalogue of PNe (1143 objects; Acker et al. 
1994), which happens to be a bulge PN. IPHAS has demonstrated 
that extremely young compact PNe can be reached (Viironen et al. 
2009b), while Sabin et al. (in preparation) have found some 20 new 
PNe with diameters of 1-3 arcsec in by-eye searches of IPHAS im- 
age mosaics. Even smaller, but brighter, nebulae around symbiotic 
stars of the dusty D subtype are emerging — the record so far be- 
ing IPHAS J193943.36+262933.1, a new D symbiotic star with an 
Ha extent of only 0.12 arcsec that has been confirmed via Hubble 
Space Telescope imaging and recently studied with the 10.4 m Gran 
Telescopio Canarias (Rodriguez-Flores et al. 2014). 

Apart from opening up new discoveries, a further benefit of good 
seeing is the clearer view of nebular structure that it offers. This 
is nicely demonstrated in Fig. 21. SHS and VPHAS-+ detect the 
main features of the PNe NGC 2438 and NGC 2899 to very similar 
depth — for example, the fainter outer halo is just detected in both 
versions of NGC 2438. But, evidently, the VPHAS+ images better 
resolve the fine sculpting within both nebulae as a consequence of 
the seeing FWHM being under a half that prevailing in SHS data. 
The extended dynamic range of VPHAS+ helps in this respect, 
too, in that early saturation also obliterates detail. This advantage is 
especially clear in the images of NGC 2899, where the structure in 
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Figure 22. A cut-out at full resolution from M8, the Lagoon Nebula. This 
is an RGB image centred on RA 18" 09™365, Dec. —24° 01’ 51” (J2000), 
and spanning 150x 150 arcsec”. The filters are combined such that R/G/B = 
Hø/i/r. 


the bright nebulous lobes is preserved in VPHAS+, but is entirely 
bleached out in SHS. The more the level of detail that can be picked 
out, the more certain and subtle morphological classifications and 
interpretations can become. 

The combination of good seeing and high dynamic range also 
makes the separation of fainter stars from background nebulosity 
much easier. This capability is critically important to the study of the 
young massive clusters, still swathed in diffuse H 1 emission, where 
the analysis of stellar content is very much a focus of continuing 
research. For example, Feigelson et al. (2013) have offered a critique 
of the nuisance created by spatially complex nebulosity. The obvious 
answer to this and the problem of dust obscuration is to turn to 
selection using NIR and X-ray data. Nevertheless, the availability 
of imaging data of the high quality seen in VPHAS+ data will 
make it possible to extend SEDs for many more stars into the 
effective-temperature (and reddening) sensitive optical domain. In 
addition, understanding the shaping of the interstellar medium in 
star-forming environments remains an important part of the picture 
(see e.g. Wright et al. 2012 on proplyd-like structures in Cyg OB2). 
The detail that the VST is capable of revealing both in obscuration 
and ionized hydrogen in star-forming regions can be quite exquisite. 
Here, in Fig. 22, we illustrate this with an excerpt from VPHAS+ 
data on the Lagoon Nebula, showing the fine tracing of the shapes 
of dark globules and eroding dusty structures that is achieved. 

In planetary and other evolved-star nebulae, it is of course im- 
portant to identify the ionizing object. The search for missing PN 
central stars is a quest that VPHAS-+ can aid greatly through the 
provision of spatially well-resolved u and g data. Indeed, inspection 
of the data used to construct Fig. 21 has revealed the probable cen- 
tral star of NGC 2899 for the first time. As shown in Fig. 23, there is 
very evidently a third very blue star just SW of the pair of stars that 
have, in the past, been scrutinized as possible companions to what 
is required to be an extremely hot (Tet > 250000 K), but probably 
faint central star (López et al. 1991). This blue object was detected 
on the night of 2012 December 20 at a provisional u magnitude 
of 18.79 +0.02. It fades through g (19.36 +0.02) to become unde- 
tected by the pipeline, and scarcely visible to eye inspection, in r. Its 
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Figure 23. The central star of NGC 2899 reveals itself. The top panel is 
a 1x1 arcmin? thumbnail of the centre of NGC 2899 as imaged through 
the u filter, while the bottom is the corresponding r thumbnail. The white 
right-angled bars pick out the position of an extremely blue, relatively faint 
star that is clearly present in all u (and g) exposures obtained, but is too faint 
for detection in r. 


coordinates are RA 09" 27™ 02872, Dec. —56° —06' 229 (J2000), 
just 1.7 arcsec from the more southerly of the pair of brighter stars 
examined before by López et al. (1991). Based on the g magnitude 
and an inferred V flux, we have determined the central star’s effec- 
tive temperature, via the well-established Zanstra method. Using 
the reddening and integrated Ha flux from Lopez et al. and Frew, 
kK. This is cooler than the temperature given by López et al. based 
on the ‘crossover’ method, but still extraordinarily hot for a central 
star well down the white-dwarf cooling track. 

It was one of the major science drivers for the merged VPHAS+ 
survey that u data, supported by g, would result in the detection of 
a broad range of intrinsically very blue objects — be they PN central 
stars, interacting binaries or massive OB and WR stars. An extreme 
example like NGC 2899’s central star provides the useful lesson 
that selection via the u — g, g — r colour—colour diagram would 
have failed to pick it out — because of the non-detection in r. In 
a case like this, the u, u — g colour-magnitude diagram has to be 
examined, in tandem with the appropriate images. 
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8 VPHAS+ PHOTOMETRY AS A REFERENCE 
SET FOR VARIABILITY STUDIES 


As the northern survey, IPHAS, has progressed over the decade since 
2003, there have been occasions on which it was possible to use the 
growing data base as a high-quality reference for checking transient 
reports — particularly of novae. The most spectacular IPHAS ex- 
ample of this was the nova and variable, V458 Vul (Wesson et al. 
2008; Rodriguez-Gil et al. 2010) where the eruption occurred a few 
months after obtaining Hw images revealing a pre-existing ionized 
nebula around the star. Indeed, there have been several instances 
in which photometry of the progenitor object has been extracted 
from the IPHAS data base and has been used to gain insight into 
the prior presence or absence of line emission or to set constraints 
on likely extinction (Steeghs et al. 2007; Greimel et al. 2012). Such 
opportunities will certainly arise with VPHAS+ - and be richer 
given the five filters offered. 

In the southern hemisphere, novae will be more frequent, as will 
other transient events. Furthermore responses to alerts, or the need 
to demonstrate long-term flux variations, can bring into use repeats 
of observations made necessary by initial quality-control failures. 
An example of this is provided by Vink et al. (2008) who used 
repeat IPHAS observations — taken on account of poor observ- 
ing conditions — to discuss the LBV candidacy of G79.29+0.46. 
With the increased attention being given to the reporting and ex- 
ploitation of transient objects (including the forthcoming Gaia 
alerts programme), this use of VPHAS+ will become more 
common. 


9 SUMMARY AND CONCLUDING REMARKS 


This paper has introduced and defined VPHAS-+, the VST Photo- 
metric Ha Survey of the Southern Galactic Plane and Bulge. The 
data taking, the rationale behind it, the data processing and data 
quality have all been described. The properties and limitations of 
the survey’s narrow-band Ha filter, NB-659, have been laid out 
and simulated in order to anticipate its performance. In addition, 
we have provided tables of the expected photometric colours of 
normal solar-metallicity stars to aid the interpretation of the sur- 
vey’s characteristic photometric diagrams — most are to be found in 
Appendix A, where the effect of changing the adopted reddening 
law is illustrated. The VPHAS+ Ha filter transmission is redder, 
wider and ~20 per cent higher throughput than its IPHAS counter- 
part — a difference that feeds through to noticeably different (r — 
Ha) colours for M stars. 

We have validated the photometry that is delivered by 
VST/OmegaCAM and subsequently pipelined at CASU, using test 
data taken of a field for which SDSS photometry is available. We 
find that the agreement is satisfactory, with the g-, r- and i-band 
calibrations differing by between 0.01 and 0.05 mag. However, for 
the time being, the pipeline calibration should be regarded as pro- 
visional — it will undoubtedly improve. Examples of the excellent 
imaging performance of the VST/OmegaCAM combination rela- 
tive to previous surveys have been provided, and we draw attention 
to the valuable archival role this first digital survey can fulfil in 
supporting discoveries of transient sources. 

Exploitation of the survey is now beginning. The detection of a 
compact ionized nebula around W26, the extreme M supergiant in 
Westerlund 1, has already been published (Wright et al. 2014). Ap- 
plications have been made for follow-up spectroscopy that will test 
the quality of selection of specialized object types that VPHAS+ 
photometry makes possible. Progress is also being made via direct 
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analyses of the photometry. For example, Mohr-Smith et al. (in 
preparation) are conducting a search for OB stars in the vicinity 
of the massive cluster Westerlund 2, and they are finding a close 
match between the properties of known cluster O stars as derived 
from VPHAS+ data and those inferred by Vargas Alvarez et al. 
(2013, see also Drew et al. 2013). This and other early appraisals 
of the data indicate that VPHAS-+ will be an excellent vehicle for 
automated searches for reddened early-type stars. Kalari et al. (in 
preparation) are employing both narrow-band Ha and the broad- 
bands to measure mass-accretion rates in pre-main-sequence stars: 
they are finding that Ha mass-accretion rates in T Tauri stars com- 
pare favourably to rates determined from the u band in the case of 
the Lagoon Nebula, NGC 6530. 

As the calibration of the survey data improves, the measurement 
of accurate integrated Ha fluxes for many faint southern PNe and 
other extended objects becomes possible, and will extend the work 
of Frew et al. (2013, Frew et al. 2014). In due course these fluxes 
can be compared with existing and also new radio continuum fluxes 
coming on stream (see e.g. Norris et al. 2011) in order to determine 
reliable extinction values for many faint nebulae currently lacking 
data. This technique has already been applied to the case of W26 in 
Westerlund 1 (Wright et al. 2014). 

When it becomes possible to cross-match VVV and VPHAS+ 
data, it will open up the power of homogeneous photometric map- 
ping of the central parts of the Galactic plane in up to 10 photometric 
bands spanning both the optical and the NIR. Beyond the VVV sky 
area, there is a synergy to be exploited in bringing VPHAS+ data 
together with those of the all-sky 2MASS survey (Skrutskie et al. 
2006) and with the UKIDSS Galactic Plane Survey (Lucas et al. 
2008), in those parts of the first and third Galactic quadrants the 
latter has covered. It is worth noting, however, that 2MASS alone is 
too shallow to link effectively with VPHAS-+ for sightlines where 
the integrated visual extinction is less than ~5 mag. This does mean 
that the longitude range 230° < £ < 300°, in particular, is presently 
lacking sufficiently deep NIR photometry. In the longer term, many 
of the sources of interest that VPHAS-+ finds will benefit from ac- 
curate parallaxes and other data from ESA’s Gaia mission — given 
the similar sensitivity limits reached. Conversely in the meantime, 
VPHAS-+ has already begun to assist ambitious wide-field spec- 
troscopy programmes such as the Gaia-ESO Survey (Gilmore et al. 
2012) through the provision of the wide-field photometry needed 
for target selection and field setup. 

By the end of 2013, 25 percent of all observations making up 
the survey had been obtained to the required quality, and in 2013 
May a first release of single-band catalogues was made to the ESO 
archive that contained roughly 10 percent of the eventual total 
(based on data obtained prior to 2012 October 15). By design, the 
characteristics of VPHAS+ are similar to those of the IPHAS and 
UVEX Galactic plane survey pair in the north. In particular, the 
double-pass strategy is shared, with the result that the majority of 
detected objects are picked up and measured twice, with no more 
than ~0.1 percent of objects missed altogether. This feature has 
informed the way in which the IPHAS DR2 catalogue (Barentsen 
et al., in preparation) has been constructed — and it is intended that 
a first band-merged VPHAS-+ catalogue, for public release, will be 
built along analogous lines during the second half of 2014. This 
will incorporate data from the first three seasons of VST observing, 
and give a complete photometric account of the Galactic mid-plane. 
For ease of use, for every detected source, the catalogue will pro- 
vide a single recommended set of magnitudes in up to five optical 
bands. 
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APPENDIX A: SYNTHETIC COLOUR REDDENING TABLES 


Synthetic colours for main-sequence stars, computed as described in Section 4, are tabulated in full in an online supplement for three 
representative reddening laws (Ry = 2.5, 3.1 and 3.8) and a range of reddenings (Ao = 0, 2, 4, 6, 8, 10). The form of the reddening laws 
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Table Al. VST/OmegaCAM synthetic colours for B, A main-sequence stars in the (u — g), (g — r) plane reddened 
with an Ry = 3.1 extinction law. (Full table online.) 


Spectral Ao = 0 Ao = 2 Ag=4 Ao = 6 Ao = 8 
type (u—g) (g-r) (u-g) (g-r) (-g) (g-r) WU-g) (g-r (u-8) (e-r) 
BOV 1.433 0.271 0.692 0.529 0.087 1.301 0.891 2.050 1.632 2.777 
BIV 1.324 0.240 0.584 0.558 0.195 1.329 0.995 2.076 1.719 2.802 
B2V 1.209 0.218 0.470 0.579 0.307 1.350 1.104 2.096 1.808 2.821 
B3V 1.053 0.186 0.315 0.610 0.460 1.379 1.250 2.125 1.923 2.849 
B5V 0.828 0.139 0.092 0.655 0.680 1.423 1.460 2.166 2.080 2.890 
B6V —0.728 —0.121 0.007 0.672 0.776 1.439 1.550 2.182 2.144 2.905 
B7V —0.580 —0.100 0.152 0.692 0.918 1.458 1.682 2.200 2.234 2.922 
B8V —0.388 —0.076 0.340 0.714 1.101 1.478 1.850 2.219 2.344 2.940 
B9V —0.198 —0.046 0.528 0.742 1.285 1.504 2.019 2.244 2.445 2.964 
AOV —0.053 —0.005 0.675 0.780 1.431 1.540 2.153 2.277 2.514 2.995 
A1V —0.019 0.005 0.709 0.790 1.464 1.550 2.181 2.287 2.525 3.005 


A2V 0.021 0.025 0.749 0.809 1.505 1.568 2.217 2.304 2.538 3.021 
A3V 0.038 0.059 0.771 0.840 1.531 1.597 2.241 2.332 2.541 3.048 
ASV 0.067 0.125 0.805 0.904 1.567 1.658 2.269 2.390 2:323 3.105 
ATV 0.044 0.199 0.788 0.975 1.554 1.726 2.252 2.456 2.474 3.169 


Table A2. VST/OmegaCAM synthetic colours for B, A main-sequence stars in the (r — i), (r — Ha) plane reddened with an Ry = 3.1 extinction law. 
(Full table online.) 


Spectral Ap = 0 Ag =2 Ag =4 Ao = 6 Ao = 8 Ao = 10 
Type (r — i) (r—Ha) (r-i) (r—Ha) (r-i) (r—Ha) (r-i (r—-He) (r-i (r—-Ha) (r-i) (r— Ha) 


BOV —0.150 0.054 0.278 0.198 0.694 0.316 1.100 0.409 1.496 0.478 1.884 0.526 
BIV —0.136 0.048 0.291 0.192 0.708 0.310 1.114 0.403 1.510 0.472 1.898 0.519 
B2V —0.123 0.045 0.304 0.188 0.721 0.306 1.126 0.398 1.523 0.466 1.911 0.513 
B3V —0.104 0.044 0.323 0.186 0.740 0.303 1.145 0.394 1.541 0.462 1.929 0.508 
B5V —0.077 0.039 0.349 0.180 0.765 0.295 1.170 0.386 1.566 0.452 1.954 0.497 
B6V —0.068 0.036 0.358 0.177 0.774 0.291 1.179 0.381 15/5 0.448 1.963 0.492 
B7V —0.057 0.029 0.369 0.170 0.785 0.284 1.190 0.374 1.586 0.440 1.973 0.484 
B8V —0.045 0.018 0.382 0.158 0.797 0.272 1.202 0.362 1.598 0.427 1.985 0.471 
BOV —0.028 0.006 0.398 0.145 0.813 0.259 1.218 0.348 1.614 0.413 2.001 0.456 
AOV —0.009 —0.005 0.418 0.133 0.833 0.246 1.238 0.334 1.633 0.399 2.020 0.441 
AIV —0.003 —0.003 0.423 0.135 0.838 0.248 1.243 0.335 1.638 0.399 2.025 0.442 
A2V 0.006 —0.004 0.432 0.134 0.847 0.247 1.251 0.334 1.646 0.397 2.033 0.439 
A3V 0.021 —0.008 0.446 0.130 0.861 0.241 1.265 0.328 1.660 0.391 2.047 0.432 


ASV 0.051 0.005 0.476 0.141 0.890 0.250 1.293 0.335 1.687 0.396 2.073 0.436 
ATV 0.083 0.027 0.507 0.160 0.920 0.268 1,322 0.350 1.716 0.410 2.101 0.448 


used is due to Fitzpatrick & Massa (2007). As an example of the tables available, we include excerpts from the second and fifth tables that, 
respectively, provide Ry = 3.1 blue-filter and red-filter colours for B, A main-sequence stars. 

Two further tables of synthetic colours are included in the supplement for K-M giants that have been computed using P98 library spectra. 
Data are provided for the R = 3.1 mean Galactic law only, for the limited purposes of (a) giving an impression of how these luminous red 
objects may contaminate the (u — g, g — r) diagrams at redder (g — r) through u red leak and (b) enabling comparisons with the M-giant spur 
commonly seen in the (r — Ha, r — i) colour—colour diagrams. 
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